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Between many prominent contributions of Markus Biittiker to mesoscopic physics, the scattering 
theory approach to the electron transport and noise stands out for its elegance, simplicity, universal¬ 
ity, and popularity between theorists working in this field. It offers an efficient way to theoretically 
investigate open electron systems far from equilibrium. However, this method is limited to situa¬ 
tions where interactions between electrons can be ignored, or considered perturbatively. Fortunately, 
this is the case in a broad class of metallic systems, which are commonly described by the Fermi 
liquid theory. Yet, there exist another broad class of electron systems of reduced dimensionality, 
the so-called Tomonaga-Luttinger liquids, where interactions are effectively strong and cannot be 
neglected even at low energies. Nevertheless, strong interactions can be accounted exactly using the 
bosonization technique, which utilizes the free-bosonic character of collective excitations in these 
systems. In the present work, we use this fact in order to develop the scattering theory approach to 
the bosonization of open quasi-one dimensional electron systems far from equilibrium. 

PACS numbers: 73.23.-b, 71.10.Pm, 73.63.Nm 


I. INTRODUCTION 

Non-equilibrium phenomena in condensed matter sys¬ 
tems are notoriously difficult to study theoretically, be¬ 
cause one is not able to rely on the universal relations 
for equilibrium state, such as fluctuation-dissipation re¬ 
lations, and because interactions, phase coherence, quan¬ 
tum and many-body effects manifest themselves in en¬ 
tirely different way compared to the equilibrium case, 
requiring for their description new techniques. Many the¬ 
oretical methods developed for non-equilibrium systems 
and phenomena can nowadays be found in textbooks, 
ranging from kinetic theories^ to the more sophisticated 
functional Keldysh technique^ However, in the context 
of the mesoscopic physics of small metallic systems, the 
scattering theory approach of Markus Biittiker to the 
electron transport^ has proven to be perhaps the most 
efficient and widely used method. The great success of 
this method may be explained by the simplicity of the 
required calculations, and by relative universality of its 
applications. Indeed, this approach relies on the free- 
fermionic character of the electron transport, justified by 
the Fermi-liquid theory of metals^ according to which 
interactions can be often neglected at low energies. 

However, the recent progress in the experimental tech¬ 
niques has revealed the new very interesting class of 
mesoscopic systems, quasi-one dimensional (ID) conduc¬ 
tors, which being perhaps less broad, become neverthe¬ 
less more and more experimentally accessible nowadays. 
Since the earlier works of Tomonaga£ and Luttinger^ 
who proposed and solved the model Hamiltonians for 
quasi-ID systems (dubbed as the Tomonaga-Luttinger 
liquids), it is known that interactions of the constitut¬ 
ing fermions cannot be neglected even at low energies, 
and very often they are not perturbative, in contrast to 
the Fermi liquid case. Therefore, a special theoretical 
technique, the so-called bosonization, has been proposed 
to deal with interactions (for a review, see Refs. [8| and 


0). The bosonization technique, roughly speaking, re¬ 
places fermions with bosons of the collective excitations, 
such as charge density or current. This nonlinear trans¬ 
formation is quite complex and non-trivial. Nevertheless, 
it is relatively well understood, rigorously described, and 
widely used for equilibrium systems of a finite size, where 
it relies on imposing periodic boundary conditions on the 
fields £ 

On the other hand, in the case of open quantum sys¬ 
tems far away from equilibrium very often one faces a 
difficulty that periodic boundary conditions cannot be 
applied to the fields. This, for example, concerns the 
mesoscopic systems, which are attached to reservoirs of 
electrons (Ohmic contacts) and voltage biased in order 
to study an electron transport, i.e., the situation that 
has been considered by Markus Biittiker in the case of 
free electrons. In order to address strong interactions 
in non-equilibrium chiral quasi-ID systems, the Ref. 0 
has proposed to solve equations of motion for the bosonic 
fields with arbitrary boundary conditions in order to ex¬ 
press fermion correlators in terms of the statistics of the 
currents of free fermionsii away from the scattering re¬ 
gion. This technique has been successfully used to ex¬ 
plain recent experiments with quantum Hall (QH) edge 
statesiiS However, the free-bosonic character of excita¬ 
tions in quasi-ID systems of fermions with a linearized 
spectrum suggests that perhaps as simple and powerful 
scattering theory as the one of Markus Biittiker can also 
be formulated in the case of strong interactions. Our 
work presents an effort in this direction. 

Earlier versions of the scattering theory for bosons in 
quasi-ID electron systems have been proposed in var¬ 
ious context, including inter-edge interactions in QH 
systems^ the universality of the DC conductance of 
quantum wires^ thermal transport^ the frequency- 
dependent linear response,— resonant dephasing in elec¬ 
tronic interferometers^ energy exchange at the QH 
edge^ equilibration of QH edge states by an Ohmic 


2 


contact^ and the decoherence of single-particle exci¬ 
tations at the QH edge^S Here, instead of focusing on 
particular physical phenomena, we formulate the scat¬ 
tering theory approach to the bosonization on a more 
rigorous level. The goal of this approach is to overcome 
limitations of the scattering theory for free fermions by 
accounting for a broad class of strong density-density in¬ 
teractions non-perturbatively. The trade-off of this tech¬ 
nique is that the fermion “mixing”, i.e., electron tunnel¬ 
ing and backscattering effects have to be taken into ac¬ 
count perturbatively^I (with an exception of the bound¬ 
ary conditions considered in Sec. El and of the electron 
mixing in reservoirs discussed in Sec. ED- This requires 
the knowledge of electron correlation functions. Our 
present work proposes a framework for the calculation 
of such correlators for open quasi-ID electronic systems. 

The rest of the paper is organized as follows. We 
start in Sec. ED with the pedagogical introduction to the 
bosonization, introduce interactions and formulate the 
scattering problem for bosons. This is followed in Sec. 
IIIII by the proof of the orthogonality and completeness of 
the basis of scattering states. This step is used to prove 
the fermionic commutation relations for vertex operators, 
and thus completes the bosonization of the interacting 
fermions. Zero modes in open systems acquire a new 
physical meaning and properties, which are the subject 
of Sec. lIVI We proceed in Sec. El with arbitrary boundary 
conditions for the fields, and connect fermionic correla¬ 
tors to the full counting statistics (FCS) of free-fermionic 
currents Jd. thereby generalizing the results of Ref. 0- 
Finally, in Sec. lVIl we combine the scattering theory with 
the quantum Langevin equations in order to account for 
the effects of dissipation and fluctuations, arising in the 
electrical circuit, to which a quasi-ID system is attached. 


II. INTRODUCTION TO BOSONIZATION 
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FIG. 1: Schematic representation of a quasi-ID mesoscopic 
system, consisting of a set of chiral fermions, interacting in 
a finite region confined inside the dotted box. Although, the 
spatial location and orientation of the channels may be arbi¬ 
trary, we choose, without loss of the generality of our anal¬ 
ysis, the same parametrization with the coordinate x. Each 
channel may originate from its own reservoir, biased with the 
electro-chemical potential p n , and having its own tempera¬ 
ture T n , n = 1,..., N. 

tation of each channel may be arbitrary. It is only for 
the convenience, and without loss of the generality, we 
parametrize the channels with the same coordinate x , as 
shown in Fig. [I] The fermions carry the charge density: 

Pn(x) = ipt(x)ip n ( x). (1) 

Using the standard commutation relations for fermions 

{i) n (x),^\,(x')} = 5 nn '8{x-x '), (2) 

one obtains the relation of the charge locality 

[pn(x),ipl,(x')\ = 5 nn ’6{x - x')ipl(x), (3) 

which imply that the operator ijr n , (a/), representing an 
electron, creates the charge equal to 1 in the nth channel 
at the point x = x'. 

The Hamiltonian of free chiral fermions with the spec¬ 
trum linearized in the vicinity of the Fermi level reads 


We start with the simple example of free chiral 
fermions, representing electrons in a quasi-ID mesoscopic 
system, and add interactions below in this section. The 
details of the bosonization procedure may be found in 
a number of textbooks, for example, in Ref. Q. The 
formulation of the problem in the context of chiral quan¬ 
tum Hall edge states, both at integer and fractional fill¬ 
ing factors, may be found, for instance, in Refs. J22| and 
[231 ]. In this section we outline, in a pedagogical manner, 
only essential steps needed for understanding the rest of 
the paper. Throughout the paper we use unites, where 

e = h = ks = 1 . 


A. Free fermions 

Let us consider a system of N chiral fermions ip n (x), 
where n = originating from their own reser¬ 

voirs with the temperatures T„, biased with the electro¬ 
chemical potentials p n . The spatial location and orien- 


H 0 = -i 



dxi/jl(x)d x ip n (x), 


( 4 ) 


where v n are the group velocities of fermions at the Fermi 
level. The equations of motion dtip n = i[Ho,ip n ], which 
immediately follow from the commutation relations (f2j), 

{d t + v n d x )%p n {x, t) = 0, (5) 


describe chiral waves that propagate with constant 
speeds. 

In equilibrium, one imposes periodic boundary condi¬ 
tions, ipn{xM) = tpn{x + L n ,t ), where L n is the size of 
the nth channel, and presents the solution as a sum over 
plane waves. Taking the thermodynamic limit, L n —> oo, 
the solution reads 


ip n {x,t) 


1 

V 2 TTOn 


OO 



—oo 


T = t-x/Vn, (6) 
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where 

{c n (u),c\,(u')} = S nn >S(uj - a/). (7) 

It is easy to see, that the operators © satisfy the com¬ 
mutation relations ©. 

Finally, we note that at ground state (c n (w)cj((u/)) = 
S(lo — — fi n ), he., all the states below the Fermi 

level are occupied, while all the states above the Fermi 
level are empty. Therefore, the free-fermionic correlation 
function reads 

We will use this result below as a reference. 


B. Free bosons 


Under the same assumption, that the spectrum of 
free fermions can be linearized in the vicinity of the 
Fermi level, namely, if the perturbations that drive the 
system away from the equilibrium are relatively weak, 
p n ,T n < £f, where £f is the Fermi energy, one may 
consider such perturbations as incompressible deforma¬ 
tions of the Fermi sea (shown schematically in Fig. [2j. 
The simplest example is given by free fermions with the 
density accumulated as a result of the shift of the electro¬ 
chemical potential: 

( Pn ) — D n fl n , E n — 1/(^7 TV n ), ( 9 ) 


where D n is the density of states of the nth channel at 
Fermi level. This relation leads to the well-known univer¬ 
sal expression for the ID charge current, ( j n ) = v n (p n ) = 

Pn/2 f . 

In the next step, we assume that these relations also 
hold locally for coordinate-dependent deformations, and 
expand the grand potential of a quasi-ID system of free 
fermions to second order in small deformations of the 
Fermi sea, 



PnQ 


n, 


( 10 ) 


where Q n = f dxp n {x) is the total charge in the nth 
channel. Varying densities, SE/5p n = 0, while keeping 
potentials p n constant, one obtains the relation © for 
the ground state. On the other hand, according to Eq. 
©, all fermions move with the constant speeds. There¬ 
fore, the densities should satisfy same equations: 


(d t + v n d x )p n {x,t) = 0. (11) 


However, if the densities are replaced by operators, and E 
is regarded as a Hamiltonian, the same equations should 
follow from the equation of motion dtp n = i[E , p n ]- This 
is possible only if the densities satisfy the following com¬ 
mutation relations: 


[p n (x), p n '(x')\ = -(i/2n)6 nn >d x 5(x - x'). (12) 



FIG. 2: Deformation of the Fermi sea in one of the chan¬ 
nels is schematically shown. All the states below the wavy 
line are filled, while all the states above this line are empty. 
This deformation leads to the accumulation of the coordinate- 
dependent ID charge density p n (x). 

This completes the quantization of the deformations of 
the Fermi sea. 

It is convenient, and quite common, to replace the den¬ 
sities and currents with the displacement fields (f> n : 

Pnipc) — 71 ^x^nip^), Jn{%) — 71 ^t^nip^), ( 13 ) 

Z7T Z7T 

where the last equation follows from the charge continu¬ 
ity relation. Then, the commutation relations read 

[^ n {x),c/) n '(x')\ = iir6 nn 'Sgn(x - x'). (14) 

We proceed, as for free fermions, by solving equations of 
motion in terms of plane waves, imposing periodic bound¬ 
ary conditions, and taking the thermodynamic limit. The 
resulting spectral decomposition for the new fields reads 

OO 

(j)n(x,t) = (fn(r) + j [e~ lUT a n {u) + h.c.] , (15) 

o 

where the creation and annihilation operators satisfy 
bosonic commutation relations 

[a n (u),aj,,(w')] = 5 nn '6(ui - a/), (16) 

and the zero modes, 

<Pu(t) = -<^ 0) - 2irv n Q n T/L n , (17) 

account for the homogeneous part of the charge density, 
and change the number of fermions Q n in each channel 
by 1, i.e., [Q„,eVi 0> ] = 1. 

Formally, the zero modes ensure correct commutation 
relations for the fields (f> n , as well as periodic bound¬ 
ary conditions. However, in the thermodynamic limit, 
L n —> oo, and for an open system, their quantum na¬ 
ture may typically be ignored (see, however, the discus¬ 
sion in Secs. m and El. Indeed, at any finite distances 
As, the quantum fluctuation of the zero mode scales as 
2TrAQ n Ax/L n ~ Ax/L n —> 0. Therefore, using Eq. © 
one can replace zero modes with their average values: 

-t {(Pn(T)) = -/i„T, (18) 

where, we recall, r = t — x/v n . For instance, according 
to Eq. (fT51) . the average stationary currents acquire the 
values (j n ) = — dtip u (t)/2f = p, n /2ir, in full agreement 
with the Fermi liquid theory. 
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C. Bosonization of fermions 

The two alternative descriptions of the chiral 
fermions outlined above are unified by the bosoniza¬ 
tion procedure^ according to which the fermions are ex¬ 
pressed in terms of bosons as follows 

ip n (x) oc e l ^ n(x) , n=l,...,N, (19) 

where the prefactor constants are determined by the 
high-energy cut-off and can be found in the end of cal¬ 
culations by, e.g., comparing the correlators of so defined 
fermions to the correlation functions (0 ■ Using the spec¬ 
tral decomposition (fTSl) . one can verify that the opera¬ 
tors (fT9l) satisfy the fermionic commutation relations 0 
and the charge locality relation 0— i The subtle step 
of connecting the densities © and ED. as well as the 
Hamiltonians (HOD and m , relies on the point-splitting 
procedure,— which accurately subtracts the contribution 
of the Fermi sea. The recent work [25| points to some 
physical aspects of this procedure in the context of the 
evaluation of the electron correlation functions and ar¬ 
rives at the conclusion that for the contributions close to 
the Fermi level this procedure can be ignored. 

Finally, we note that the correlation functions 0 of 
free fermions in the ground state may now be obtained 
by using the operators (fTflll , substituting the spectral de¬ 
composition of the bosonic fields ED, and applying the 
normal ordering procedure. Importantly, the commuta¬ 
tion relation for zero modes, [Qn,e lv ’» ) ] = 1, gaurantees 
the independence of fermions that belong to different ID 
systems of finite size L n : (i/> n i/d,) = 0 for n ^ n'. This 
follows from the orthogonality of quantum states of fi¬ 
nite ID systems containing different number of fermions. 
However, for open quasi-ID fermionic systems connected 
to reservoirs this argument has to be reformulated, which 
is discussed in the Sec. eh 

D. Interactions 

Rewriting the total energy of free fermions in Eq. m 
in terms of bosonic fields, one obtains the bare Hamilto¬ 
nian 

H 0 = ^ v n f dx \d x (j) n {x)} 2 . ( 20 ) 

77 n 

We are interested in the situation, where the interactions 
between fermions are present in the finite region of space, 
as schematically shown in Fig. [Q We wish to consider 
strong Coulomb interactions. However, in real systems 
Coulomb interactions are screened in a quite complex 
way. Therefore, in order to keep generality of the fol¬ 
lowing analysis, we present the interaction part of the 
Hamiltonian in the general form 

H i = JI dxdyU nn >(x,y)d x <f> n (x)d y (l) n ’(y ), 

n,n' 

( 21 ) 


where the kernel U nn ' is an arbitrary potential of the 
density-density interaction, which accounts spatial sepa¬ 
ration of channels and screening effects. 

The equations of motion for the bosonic fields, = 
i[Ho+Hi, </>„], immediately follow from the commutation 
relations ED. 

(d t + v n d x )(j) n (x,t) = ^2 J dyU nn i(x, y)d y (f> n i (y,t). 

n' 

(22) 

These equations may be accompanied with the non¬ 
trivial boundary conditions, as described in Sec. El and 
solved directly. However, as a first step, we wish to ex¬ 
press fields in the second-quantized form. 

We have done so for free-fermion case, U nn > = 0, where 
the fields can be expanded in terms of plane waves, see 
Eq. (ED- Taking into account our assumption, that the 
interaction is localized in the finite region of space, we 
look for the solution of the equations (E21) in the form 

4> n (x,t) = if n (x,t) + 5(/> n (x,t). (23) 

where the zero mode contribution reads 

<Pn(x, t) = -y n t + <p n (x) (24) 

and the fluctuating part is decomposed in the oscillator 
modes as 

OO 

S(j>n(X, t) = f ^2 [^mnu(x)e~ lult a m (uj) + h.c.] . 

0 m 

(25) 

Here, the operators a m (oj) are defined in Eq. ED. and 
the wave functions &mnui(x) are the scattering states, 
which acquire the form of the plane waves far away from 
the scattering region. In the next section, we develop 
the formal scattering theory and prove the orthogonal¬ 
ity and completeness of the scattering states in order to 
guarantee the bosonic commutation relations (THl) of the 
operators ED- The zero modes (f n (x, t) deserve a special 
consideration, which is done in Sec. hyi 

Finally, assuming the Gaussian character of fluctua¬ 
tions, one can expand fermionic operators ED to second 
order in the fields 6<j> n , average over the bosonic states, 
and re-exponentiate the result. This procedure leads to 
the following expression: 

,t')} = -iy n (t - t') 

+ i[<p n (x) - (p n {x')\ + G n {x, x', t - t'). (26) 

Here, first term on the righ hand side may be attributed 
to the energy shift due to the applied chemical potential 
fi n , while the second term is the Friedel phase shift in¬ 
duced by the interaction, since d x (p n {x) = 2i r(p(x)). The 
fluctuation contribution 

G n {x, x ', t) = ([ Sc/) n (x , t) - Scj) n {x, 0 )]54> n (x', 0)) (27) 

may be found by substituting the spectral decomposition 
ED an d taking the principle value of the integral over 
co, because the co = 0 contribution has been attributed 
to the zero mode. 
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III. SCATTERING STATES 


of scattering matrix to arrive at the following expression 


In this section we formulate the scattering theory based 
on equations (l22l) . Namely, we introduce scattering states 
and prove their orthogonality and completeness. In our 
case, the scattering states are represented by N sets of 
functions, ^ mn u(x), m,n = 1,..., TV, where the first in¬ 
dex enumerates the sets, and the second index enumer¬ 
ates the functions in a particular set. They satisfy the 
equations of motion (1221) , which in the Fourier space read 


w 

I dx ^ ' ^l nU j'dx*&rnnu — kJ(5/ m 
-w 

e i(k' m -k m )W _ e -i(k' m -k m )W 


As W —> oo, the last term in this equation becomes a 
(5-function, and we arrive at the orthogonality relation: 


(iui 


V n ) ^rnnuj (^) 



dyU nn > (x,y)dy$ mn'uj (y), 


(28) 


and the following boundary conditions at \x\ -A oo, 


dx^2 = —2niu5i m 5(u' - u ). (32) 

J on 


In order to prove the completeness of scattering states, 
let us multiply Eq. (1321) by the function 4> mra / w (:c'), sum 
over m, and integrate over u. The result reads 


= [5 mn 8(-x) + S nm (uj)8(x)] e lknX , (29) 

where k n = ui/v n , and S n m are the elements of the scat¬ 
tering matrix. In other words, the mth scattering state 
describes one incoming mode in the mth channel, and 
outgoing modes in all N channels. 

As a first step, we wish to derive a useful formal re¬ 
lation for scattering states. Let us consider the sets of 
functions &mnw( x )i which satisfy Eqs. (1251) with u on 
the left hand side replaced by —u. We multiply equations 
of motion for the functions $ > *nnu( x ) by d x ^i nu] '{x), and 
equations for the functions ^i nu '(x) by d x ^^ nnu {x), inte¬ 
grate over x and sum over n, and subtract one result from 
another. Due to the symmetry of the interaction poten¬ 
tial, U nn >(x,y) = U n ' n (y,x), the right hand side of the 
equations cancels, and we arrive at the following equa¬ 
tion: f dx'E n (ud x $i nu >$t nnu + u'$i nu >d x $’Z nmil ) = 0. 
It is convenient to cut this integral at large distances, 
\x\ = W, beyond the interaction region, and integrate by 
parts. Using the asymptotic form (T2U1) . we obtain: 

w 

(uj’-uj) f dx' 52 $ lnu .d x $* rnnu =uj 6 lm e-^- k ^ w 
-w n 

-u^SniWi^uy^-^w. (30) 

n 

We note that, by choosing in this equation u/ = u, we 
immediately arrive at the unitarity of the scattering ma¬ 
te, T,n S nl S nm = 5 lm- 

Next, we extend the integral in Eq. (13U1) to infinity: 
W —^ oo. Then, for u’ ^ u the right hand side of the 
equation is fast oscillating function, which vanishes upon 
coarse graining, leading to the orthogonality of scattering 
states. On the other hand, a care has to be taken when u/ 
approaches u. In this case, we may rely on the unitarity 


OO OO 

/ duj 

u>( x ) 

0 

=-27T i$i n > u '(x'), (33) 

i.e., the expression in square brackets is the unity opera¬ 
tor in the space of chiral scattering states. Similarly, mul¬ 
tiplying Eq. (1321) by the function 4>* nn i u (x ') and repeating 
the above steps, we arrive at the analogous expression for 
the unity operator in the space of the conjugated (anti- 
chiral) states. By combining these two expressions, we 
obtain the completeness relation: 

OO 

j ^-[^ n , u ( x ')d x ^ mnu (x) 

m J Q 

~ ®mn'w(%')dx®mnu( x )] = S nn '5(x - x'). (34) 

One may easily check this relation with the simple exam¬ 
ple of U nn ' = 0, where $ mnw (a;) = 6 mn e lknX . 

To conclude this section, one of our main results is 
that, as one can easily check, the commutation relations 
0 for the fields C5l) follow from the relations of the 
completeness (13H) . One can extend the results of this 
section to the case, where M of IV available channels 
are compactified on the circles of finite size L n with the 
periodic boundary conditions for the fields. This lim¬ 
its the number of scattering states to N — M remaining 
unconfined channels. The compactified channels do not 
contribute to the right hand side of the equation (1301) . 
because the boundary terms cancel after integration by 
parts. On the other hand, on the left hand side of the 
orthogonality relation (1321) the sum over n and integral 
over x split in two parts corresponding to the two sorts 
of channels. Finally, the completeness relation (1341) does 
not change its form, however, the sum runs over scatter¬ 
ing states, i.e., from n = M + 1 to n = N. 


J dx^2 $ Inu .>' (x) 
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IV. ZERO MODES 

Let us recall, that we consider an open system, there¬ 
fore zero modes in most cases may be considered classi¬ 
cal fields, which satisfy equations of motion E?l) . Below, 
however, we will quantize zero modes in order to account 
for the quantum effects of a circuit, to which the system 
is attached. Away from the scattering region, x —> —oo, 
zero modes obviously acquire the form (1151) found ear¬ 
lier for a translationary invariant system. Therefore, we 
are looking for the solution in the form (|24| . where the 
coordinate-dependent term satisfies the equation 

Vnd x (f>n{x)+^z Y, J dyU nn >(x,y)d y & n '(y) = /z„. (35) 

n' 

This equation may be interpreted as a condition of con¬ 
stant potential at the channel n, where the first term 
on the left hand side is the contribution from the finite 
compressibility of the Fermi sea (Thomas-Fermi correc¬ 
tion), while the second term results from the Coulomb 
interactions. 

Such electrostatic problem can be formally solved with 
the help of the so called “characteristic” potential s 26 ! 27 
fmnix). Namely, one can write zero modes in the form: 

ip n {x, t) = -Unt + y ^(lXm/Vm)fmn(x), (36) 

m 

where the characteristic potentials satisfy the following 
equations 

d x fmn(x) + — J dyU nn ' (a?, y)dy fmn' (y) — (5 mn • 

n' 

(37) 

The boundary conditions may be fixed with the help of 
Eq. ifTHl) . so that the asymptotic forms read: 

fmni.x') — x8 mn T ^(x)A i /* mn , at |x| y oo, (33) 

where A f mn are the interaction-induced phase shifts. 

The equations m for the characteristic potentials 
may be solved directly, e.g., perturbatively with respect 
to the potentials U nn >. However, if scattering states 
$ mnu (j:) are already known, one can, alternatively, ex¬ 
tract characteristic potentials by evaluating the limit: 

r / \ i. *&mnu>\X) fimn /on '\ 

fmn(X) = lim ---. (39) 

/c m —>-0 % rvjYi 

This simply follows from the fact that the expression on 
the right hand side of this equation satisfies the equations 
(EH for the characteristic potentials, which can be eas¬ 
ily seen from Eqs. (1281) . Moreover, in the low-frequency 
limit and at x —> —oo we find [<f 'mnui{x) — Smn]/ikm = 
(S m „(e ,lmI — 1 )/ik m —► x6 mn , i.e., the expression (l39l) 
satisfies the boundary conditions for the characteristic 
potentials. Finally, according to the Eq. E51) . the phase 
shifts may be found from the scattering matrix: 

A fmn = lim - - -. (40) 

fc m ->0 lk m 


We stress that, in contrast to one-dimensional low-energy 
fermions, the bosons do not scatter at low frequencies, 
i.e., SV^O) = 5 mn , and the limit CUl) is well defined. 

The equations (137>1) for zero modes and © for the 
characteristic potentials may be used in order to evaluate 
voltage dependent and interaction induced specific phase 
shifts, which contribute to electron correlation functions. 
Such phase shifts are relevant for a number of experi¬ 
mental situations, as has been demonstrated in Ref. [22| . 
However, so far we have considered an ideal situation, 
where the system is attached to an electrical circuit via 
voltage biased ohmic reservoirs. Sometimes, one is not 
able to neglect effects of fluctuations inside an electrical 
circuit, which may propagate down to the system along 
the ID channels. A formal way to account for these ef¬ 
fects is to impose apropriate boundary conditions on the 
boson fields, as we demonstrate in the section ED How¬ 
ever, in the case, where the characteristic frequency of 
circuit fluctuations, oj o, is much smaller than the the in¬ 
verse time of the propagation of fluctuations through the 
system, 1/f/, there is an alternative way to proceed. 

To address this specific situation, let us consider a low- 
frequency limit of the oscillator part of the fields (f> n in Eq. 
(| 25 f) . which accounts for the dynamical effects. By using 
the equation (1551) . we arrive at the following expression 

= —6n(t) + y ^y- 1 fmn(x)d t 0m.(t), (41) 

m 

where we have introduces the operator 

OJ c 

9 n (t) = - J [e-^a n (oj) + h.c.] , (42) 

o 

and u> c is the energy cutoff, such as ujq <C w c 1 /tf. 
Note, that ( dt9 n ) = fj , n , and if the fluctuations can be 
neglected, we come back to the equation (l36l) . In order 
to take into account fluctuations, the operator dtd n has 
to be considered a circuit variable, namely, a fluctuating 
potential in reservoirs. Note also, that Eq. m repre¬ 
sents the first two terms of the adiabatic expansion. In 
particular, it is easy to see that the second term is small 
as compared to the first one by the parameter uiotf -C 1. 
However, it accounts for the interaction effects, not con¬ 
tained in the first term. 


V. BOUNDARY CONDITIONS AND FULL 
COUNTING STATISTICS 

We have already mentioned, that it might be useful in 
a number of situations to express the bosonic fields cf) n 
in terms of their values away from the scattering region, 
where the statistics of their fluctuations is assumed to 
be known. Then, by solving equations of motion (1221) 
with corresponding boundary conditions, one can find 
correlation functions of the fields 8<j) n . This method has 
been proposed in Ref. M and successfully applied to a 
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number of physical phenomena in quasi-ID systems far 
from equilibrium. 12 - Here we generalize this method to 
the case of arbitrary scattering. 

Let us assume that at the distance W upstream the 
scattering region (see Fig. [[} the fields 5<j> n are known. 
According to the asymptotic form (l29l) of the scattering 
states, and taking into account the decomposition (E5ll . 
we can write 





oo 



0 


e -M w /v*+t) an ( u ) + hc 


(43) 


On the other hand, according to Eq. m the fluctuation 
of the charge injected into the system through the cross- 
section x = — W is equal to 


L 

Sq n (t) = j dt! <5 j n (-W,t') = --^-5(j> n (-w,t). (44) 

—OO 

Comparing these two equations, one finds that 

OO 

a n (w) = -yfre luW/v " J dte iut 5q n {t ), w > 0. (45) 

— OO 

By substituting this expression into Eq. we finally 
obtain 

OO 

S(f> n {x,t) = - I dt'^2^ mn (x,t-t')Sqm{t'), (46a) 

-oo 

oo 

®mn(x,t-t')= J du^ mnu] {x)e~ lul( ' t ~ t \ (46b) 


Here, we have used the property [& m nw(x)\* = 

$ mn _ u (r) and dropped the phase factor e luW / v ™, which 
merely shifts time (note, that the correlation functions 
for stationary processes do not depend on such time 
shifts). Omitting this phase factor simply amounts to 
redefining the scattering state. 

We note that the statistics of the fluctuations of the 
fields (f> n {x, t), which are needed in order to calculate the 
fermion correlation functions, may be expressed, using 
Eqs. (USD, in terms of the statistics of fluctuations of the 
charges q n (t) transmitted through a given cross-section. 
Thus, the evaluation of fermionic correlation functions 
reduces to finding the full counting statistics of the trans¬ 
port of free fermions.— Importantly, this method allows 
one to relax the assumption of the Gaussian character of 
fluctuations needed in order to arrive at the result (26} 
and (27} . 


VI. LANGEVIN EQUATIONS 

In Secs. Hi] and mi] we have considered conservative sys¬ 
tems, while the results of Secs. m and 0 may also be 


FIG. 3: One of the reservoirs of electrons (an Ohmic con¬ 
tact) is schematically shown on the left. On one side, it is 
attached to an electrical circuit and receives the fluctuating 
current 8j c ,n from it, shown by the dashed line. This leads 
to fluctuations of the potential, Sfj, n . On the other side, it is 
connected to a quasi-ID electron system via two chiral chan¬ 
nels carrying one incoming current 5ji n ,n and one outgoing 
current 8j ou t,n- On the right, and example of a four-terminal 
systems is schematically shown. 

used to account for the dissipation in electrical circuits, 
or, e.g., in tunnel junctions— However, as far as cir¬ 
cuit effects are concerned, it has been suggested earlier 
in Ref. [19], that perhaps the most natural and efficient 
way to account for the dissipation is to apply the method 
of quantum Langevin equations. This method is based 
on the observation that electrical circuit elements typi¬ 
cally create only Gaussian fluctuations, and the system 
in the bosonic sector remains Gaussian. Therefore, when 
accounting non-linear effects (such as weak tunneling or 
weak backscattering) perturbatively, one can describe the 
dynamics of the fields by using linear equations. 

Let us consider N electron reservoirs, in general at 
different temperatures T n and different potentials /i n . 
These reservoirs are connected, on one side, via an elec¬ 
trical circuit, characterized by the frequency-dependent 
conductance matrix G nrn (u>). On the other side, the 
reservoirs are attached to a quasi-ID electron system via 
chiral electron channels. This situation is illustrated in 
the Fig. [3j where one of the reservoirs is schematically 
shown on the left. For simplicity only, and to illustrate 
our idea, let us assume, that each reservoir absorbs one 
incoming electron channel, and emits one outgoing chan¬ 
nel (in the context of the QH effect this situation corre¬ 
sponds to the case of filling factor 1). 

We first concentrate on the fluctuation contribution to 
the fields (l23l) and note that, according to the scatter¬ 
ing theory [see Eq. (HS1) ]. it is determined by the cur¬ 
rents <5jout,n outgoing from the reservoirs. Thus, one can 
start with the equation for the charge conservation in the 
Fourier space, 

itoCndUn = 8jc ,n T ^Jin,n ^'jout.n; (47) 

where C n is the charge capacitance of the nth reservoir^ 
5j c ,n is the fluctuation of the current incoming from the 
circuit, and the last two terms are the contributions of 
the electron channels. 

The currents have contributions from the fluctuations 
of the collective modes, as well as from the Langevin 
sources. The current from the circuit acquires the fol- 









lowing from: 


5jc,n = ^2 G nmSfi m + Sj a c , n , ( 48 ) 

m 

where Si* v is the source. Similarly, as it has been shown 
in Ref. [19j, the outgoing current in the electron channel 
has two contributions, 

^Jout.ra = Gq5fl n + Sj out n , (49) 

where G q = e 2 /27r/i is the conductance quantum (restor¬ 
ing physical unites), and 6j* ut n is the equilibrium ID 
current source originating from the reservoir. 

Finally, rewriting Eqs. (1461) in the frequency domain 
and in terms of currents injected to the system 

$4>n{x , w) = -(27 ri/u) ^2 ®mnu(x)6j out ,m(u) (50) 

m 

and using Eq. (1291) . one connects incoming and outgoing 
ID currents 

&jin,n(^') = ^ ' t ^nm(^)^Jout,m(^ J ); (51) 

m 

where A/i n ,n(—w) = SjJ n n (ui). These equations complete 
the set of the equations that should be solved for the cur¬ 
rents djout.n in terms of the sources. Knowing currents 
Sjont,n, one can find the fluctuating part of the fields 5<p n 
using Eq. (|50jl . However, this has to be done with cau¬ 
tion, because the matrix of conductances is degenerate. 
One can simply assume that one of the electron reservoirs 
is grounded, and corresponding potential does not fluc¬ 
tuate, or, alternatively, one may add an extra grounded 
electrode. 

Since fluctuations are assumed Gaussian, it remains to 
accompany these results with the two-point correlation 
functions of the sources. It is natural to assume, that the 
sources are at local equilibrium, and correlation functions 
satisfy fluctuation dissipation relations: 

(Sf 0Uttn (u)6f 0Ut , m (u ;')> = 5 nm 5{u+u') i (52) 

and 

(5j s c , n M<5jcV(^)) = + *') (53) 

where because of the chirality of ID electrons the factor 
of 2 in the first equation is missing as compared to the 
second one, and in the second equation we assumed that 
the circuit is at equilibrium with the bath at the tem¬ 
perature T. After the fields are expressed in terms of 
the sources, one can use Eqs. (1521) and (1571) to find the 
correlators of the fields (l27l) . and eventially, the electron 
correlation functions (l26l) . 

We conclude this section by addressing the zero mode 
contributions ip n (x,t ) to the fields cj) n (x,t). It has been 
emphasized in Sec. Ill Cl that the components of zero 


modes guarantee the independence of fermions belong¬ 
ing to different ID systems of finite size, and that in 
open quasi-ID systems this argument has to be reformu¬ 
lated. Here we propose to consider the whole system, 
i.e., the set of ID fermion channels plus an electrical cir¬ 
cuit, as an isolated ensemble of fermions, in which the 
number of particles does not change. This allows one 
to omit the phase factors ‘ in the bosonized repre¬ 
sentation of the fermionic operators, because fermions 
from different channels belong to the same system, are 
mixed (scattered) in the electrical circuit, and there¬ 
fore are not strictly independent. Nevertheless, in prac¬ 
tice their statistical independence emerges from the fact, 
that different channels are connected to each other via 
reservoirs, where the phase fluctuations are strong. For 
example, if two fermionic channels are connected via 
one reservoir, the corresponding correlator behaves as 
In(V’nV’n') ^ —Tntd, where T n is the temperature, and 
td is the dwell time of the fermion in the reservoir. Thus, 
this correlator vanishes in the thermodynamic limit. 

The remaining components of the zero modes can be 
found by using the equations (1551) and (1371) . According 
to the equation (l36l) . the zero modes are determined by 
the potentials of the reservoirs p n , which in turn simply 
follow from the Kirchhoff’s law, 

( jc,n ) T (jin.n) (jout.n) = 0. (54) 

These equations have to be accompanied by the dc ver¬ 
sions of the equations (1751) . (1791) . and m : 

(jc,n) = E (~*nm (0)/i 

mi Oout,n ) — G q p ni 
m 

(jin.n) = E Snm (9) (jout.m) (55) 

m 

and solved for p n . Here, in the notations of this section, 
'S'nm(O) = 0,1 is simply a connectivity matrix, because 
plasmons do not scatter at zero frequency. 

VII. DISCUSSION 

In the end, we would like to make important remarks 
and give practical recommendations concerning the ap¬ 
plication of the scattering theory. First of all, although 
we consider chiral systems throughout the paper, this is 
done for the convenience in order to simplify the deriva¬ 
tions. With some limitations^ on otherwise quite broad 
class of long-range interactions, our approach also ap¬ 
plies to non-chiral systems, such as Luttinger liquids con¬ 
nected to free-fermionic reservoirs. In this case, the reser¬ 
voirs are modelled by gradually switching off the inter¬ 
action at the interface with the Luttinger liquid, which 
leads to the Andreev-type process, restoring the univer¬ 
sality of the conductance of such systems*^ Therefore, 
the interaction remains localized in a region of finite size, 
while incoming and outgoing states are still free. This is 
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exactly the situation, where our method applies. Alter¬ 
natively, one can simply change the basis of scattering 
states to the left and right movers and consider the in¬ 
terfaces with reservoirs as additional scattering centers. 

Second, in QH systems at integer filling factors larger 
than 1, co-propagating electron channels interact even in 
the asymptotically remote regions. However, this inter¬ 
action is homogeneous and can be easily diagonalized by 
applying rotations of the basis in each sector of incom¬ 
ing states belonging to the same edge. Thus, the bosonic 
fields appearing in the vertex operators (fTi?l) can still be 
expressed in terms of actual scattering states. 

Third, additional electrostatic potentials, V n (x), e.g., 
induced locally by metallic gates, can easily be accounted 
by adding linear terms (l/27r) / dxV n d x <j) n to the 

Hamiltonian Hq + H±. They slightly modify the equa¬ 
tions of motion (l22l) by adding the term —V n (x) to the 
right hand side. Consequently, this term arises in the 
right hand side of Eq. (1551) for zero modes, modifying the 
corresponding electrostatic problem, and shifting densi¬ 


ties: (1/27 T)d x (f> n —>■ (l/27r)d x tf>„ + A p n . In turn, this in¬ 
troduces additional phase shift 2n f dxAp n in fermionic 
operators m, in agreement with the Friedel sum rule. 

Finally, we note that it would be interesting to ex¬ 
tend our scattering theory in the analogy to the Floquet 
theor y 29 ! 30 in order to investigate the photon-assisted 
electron transport in quasi-ID electron systems. This 
seems to be an obvious and straightforward next step, 
because even if a system is biased with time-dependent 
potentials, it remains Gaussian, and therefore the bosons 
are free. 
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